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15 Abstract 

Purpose: Photoacoustic imaging has proven to be able to detect vascularization-driven optical 
absorption contrast associated with tumors. In order to detect breast tumors located a few cen- 
timeter deep in tissue, a sensitive ultrasound detector is of crucial importance for photoacoustic 
mammography. Further, because the expected photoacoustic frequency bandwidth (a few MHz to 

20 tens of kHz) is inversely proportional to the dimensions of light absorbing structures (0.5 to 10+ 
mm), proper choices of materials and their geometries and proper considerations in design have to 
be made to implement optimal photoacoustic detectors. In this study, we design and evaluate a 
specialized ultrasound detector for photoacoustic mammography. 

Methods: Based on the required detector sensitivity and its frequency response, a selection of 

25 active material and matching layers and their geometries is made leading to a functional detector 
models. By iteration between simulation of detector performances, fabrication and experimental 
characterization of functional models an optimized implementation is made and evaluated. For 
computer simulation, we use ID Krimholtz-Leedom-Matthaei (KLM) and 3D finite-element (FEM) 
based models. 

30 Results: The experimental results of the designed first and second functional detectors matched 
with the simulations. In subsequent bare piezoelectric samples the effect of lateral resonances was 
addressed and their influence minimized by sub-dicing the samples. Consequently, using simu- 
lations, a final optimized detector was designed, with a center frequency of 1 MHz and a -6 dB 
bandwidth of 0.4-1.25 MHz (fractional bandwidth of ~80%). The minimum detectable pressure 

35 was measured to be 0.5 Pa. 

Conclusion: A single-element, large-aperture, sensitive and broadband detector is designed and 
developed for photoacoustic tomography of the breast. The detector should be capable of detecting 
vascularized tumors with 1-2 mm resolution. The minimum detectable pressure is 0.5 Pa, which 
will facilitate deeper imaging compared to the current systems. Further improvements by proper 

40 electrical grounding and shielding and implementation of this design into an arrayed detector will 
pave the way for clinical applications of photoacoustic mammography. 
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I. INTRODUCTION 



Photoacoustic imaging is an intrinsically hybrid biomedical imaging modality which is 
based on light excitation and ultrasound detection. In photoacoustics, short-pulsed laser 
light diffusively penetrates into tissue, is selectively absorbed by specific chromophores, 
such as hemoglobin, causing ultrasound waves to be generated by thermoelastic expansion. 
These are then detected by ultrasound detectors at the surface of tissue. Since ultrasound 
scattering is considerably lower than light scattering in biological tissue, the information 
carried by ultrasound waves arrives with less losses and distortion at the tissue surface than 
light waves would do. Thus photoacoustic imaging takes the advantages of both optical 
imaging and ultrasound imaging, with a high optical contrast, a high ultrasound resolution 
and large imaging depths |l|-l6|. One of the most important applications of photoacoustics 
is breast imaging, where the optical contrast comes from higher average hemoglobin levels 
associated with malignant masses compared to healthy breast tissue [71, |8|. The technique 



promises to be an a. 
breast cancer 



14. 



ternative to ionizing x-ray and low contrast ultrasound imaging to detect 



The faithful and sensitive detection of ultrasound (US) lies at the heart of a photoacoustic 
imaging system. The US detector largely determines image contrast, resolution and imaging 
depth of the system. In photoacoustic (PA) imaging of the breast, there is a requirement 
to detect tumors located a few centimeter deep in tissue, where hght is heavily attenuated. 
Thus a sensitive ultrasound detector is of crucial importance. Further, the frequency ran ge o f 
photoacoustic waves is inversely proportional to the dimensions of absorbing structures 15 |. 



In breast tissue, structures of interest range from 0.5 to 10-|- mm initiating ultrasound 
frequencies from a few MHz down to tens of kHz {l^. Thus a broadband US detector is 
required, centered on an optimum frequency. 

Various US detectors with different specifications have been employed in photoacoustic 



(thermoacoustic) systems for breast imaging 10 



-12 



15 



17|-|20|. In 2000, Kruger et al. [17|, 



reported the first thermoacoustic (TA) breast scanner in CT configuration. The US detector 
used 64 commercially available immersion detectors (1-3 piezocomposite, model 3847, Pana- 



metrics, Waltham, Mass.) each with flat face and a diameter of 13 mm, arrayed in a spiral 
pattern spanning the surface of a hemispherical bowl. The detector had a center frequency 
of 1 MHz with a fractional frequency bandwidth (FEW) of 70%. The sensitivity of the de- 
tector was not reported. A later version of the system from 2010 jlo|, employed an array of 
128 detectors (1-3 piezocomposite), 3 mm in diameter with a 5 MHz center frequency and a 
FEW of 70%. Submillimeter breast vasculature down to a depth of 40 mm was successfully 



visualized 



10| . The sensitivity was not reported. However, the authors recommended the 



use of detectors with a lower center frequency to have greater imaging sensitivity. 

n 

Pramanik et al. [12] in 2008, reported a breast cancer detection system combining TA 
and PA tomography. They used 13-mm/6-mm-diameter active area non-focused detectors 
operating at 2.25 MHz center frequency (piezocomposite, ISS 2.25 x 0.5 COM, Krautkramer) 
for signal detection. The detectors used have a large active surface area and a relatively low 
center frequency to gain a high sensitivity. The FEW is reported as varing between 60-120% 
but no details of sensitivity are reported. 

The laser optoacoustic imagin g sy stem (LOIS) was developed over a decade ago and has 

uses an array of 64 wideband 



undergone several iterations 
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polyvinylidene fluoride (PVDF) (bandwidth upto 2.5 MHz) elements arranged in a concave 
arc. Each detector element has a rectangular surface, with a large aspect ratio (20 mm x 3 
mm). This design gives the detector a slice-shape focusing area, providing high sensitivity in 
the region of interest. The sensitivity of this system was reported to be 1.66 mV/Pa at 1.5 
MHz, the minimum detectable pressure (MDP) was not reported ll|. However, for some 
earlier versions the MDP values have been presented including measured {21 1 and estimated 
using only the thermal noise generated by the detector capacitance according to the Nyquist 
law [id, 122! . 

The Twente Photoacoustic Mammoscope (PAM) has been previously developed in our 



group in 2004 



24j . US detection uses a planar array of 590 PVDF detector elements. Each 



element has a 2 mm x 2 mm active surface area. The detector has a center frequency of 1 MHz 

n 

and a fractional bandwidth of 130% [20|. The measured MPD value is 80 Pa [9|]. Promising 



I4J . |25|. However, the 



clinical measurement results have been reported in 2007 and 2012 
system still suffers from the relatively low sensitivity of the PVDF detector. 

PVDF detector arrays give LOIS and PAM broad bandwidths, while the measured sen- 
sitivity of the detectors are generally low compared to PZT detectors due to the reduced 



electromechanical coupling coefficients [235. However, PZT detector suffers from a relatively 
105 low bandwidth. As a compromise, piezocomposite detectors provide better sensitivity and 
reasonably good FBW as employed by Kruger et al. 10(] and Pramanik et al. 12|]. 

To achieve a high sensitivity of the detector, we choose a highly sensitive PZT material, 
while tailoring the bandwidth of the detector to be reasonably wide. A single element PZT 
detector structurally consists of the active piezoelectric material, front- and back-matching 
no layers and a backing layer. To have both high sensitivity and broad bandwidth, the materials, 
their acoustic characteristics and their dimensions should be carefully chosen. Furthermore, 
the aperture of the detector should be o pti mized as there is a trade-off between the lateral 
resolution and sensitivity of the system [26 1. 

In this paper, we present the design considerations for this specialized ultrasound detector 
115 for PA breast imaging. We specify the most important detector output characteristics such 
as sensitivity and frequency response, and justify the selection of active material and match- 
ing layers and their geometries. We iterate between simulation of detector performance, 
fabrication and experimental characterization of functional models to arrive at an optimized 
implementation. For computer simulation, we use ID Krimholtz-Leedom-Matthaei (KLM) 
120 models and 3D finite-element (FEM) based models. 



II. DESIGN PARAMETERS 

A. Sensitivity and acceptance angle 

In photoacoustic tomography, the breast is illuminated with short-pulsed laser light and 
the ultrasound detector scans the object ideally through 360°. A requirement in this geom- 

125 etry is that each US detectors's acceptance angle is wide enough to detect photoacoustic 
signals generated throughout the entire object for each detector angular position around 
the object (Figured]). During a backprojection style reconstruction, the coherent signals 
from all detector angular positions are summed up to form an image of the photoacoustic 
sources [27^. For a single element square surface detector, the -6 dB acceptance angle (0) 

130 depends on the width {W) and center frequency (/o) of the detector and can be expressed 

5 



as 



28| 



7Q 

^ = 2sin-(^) (1) 

with k the wavenumber. Equation 1 shows that the acceptance angle decreases as the 
aperture increases. However a larger surface area is important since for the same piezo 

135 element thickness, a lower MDP is attained. This is due to a lower thermal-induced noise due 
to the higher electrical capacitance associated with larger detectors. To use a sensitive large 
aperture detector without compromising the acceptance angle, Li et al [291] and Pramanik 
et al [30] introduced the concept of a negative acoustic lens. Here an acoustic lens enlarges 
the acceptance angle of a large aperture detector improving the lateral resolution of the PA 

140 tomographic system. This provides a solution where one can use a large area detector while 
maintaining a wide acceptance angle. 

Based on the above considerations, the choice is made for a large aperture detector with 
a square surface (5 mm x 5 mm) with an appropriate lens 3l|, which will be discussed in 
a future article. The square shape is chosen for convenience in array development using 

145 traditional dicing. 



B. Center frequency 



Photoacoustic breast imaging is based on the tissue optical contrast due to vasculariza- 
tion. This is enhanced around a tumor due to the process of angiogenesis 32|. This process 
is reported to go through two phases separated by the "angiogenic switch". Exponential 
150 tumor growth ensues in the second phase (vascular phase), which occurs from tumor sizes 



33| . This indicates that the resolution of our system is preferably to 



of 1-2 mm in diameter 
be smaller than 2 mm. 

The upper and lower limit of the frequency range of the detector {fmax and /m.m(MHz)) 
can be estimated knowing the smallest and largest sphere {amin and amax (mm in diameter)) 



155 that are required to be resolved by the system using 



19|: 



'■^ J rain 



(2) 



in which v (mm /is ^) is the acoustic velocity of the medium. Therefore, faithful registration 
of spherical tumors with diameters from 2 mm to 10 mm requires an ultrasound detector with 



frequency bandwidth from 2.25 MHz down to 48 kHz. To achieve the required bandwidth, 
160 the center frequency of the detector is designed to be 1 MHz, while optimizing the bandwidth 
of the detector to be reasonably broad. 



III. MATERIALS AND FABRICATED MODELS 



A. Materials 



As mentioned earlier, we prefer PZT detectors due to the superior dielectric constant. 



165 lower dielectric loss and higher coupling coefficients j23| . A high sensitivity commercial piezo- 
electrical material CTS 3203HD (CTS Communications Components, Inc., Albuquerque, 
NM) is used for the active material due to the higher coupling coefficients compared to 



other common used commercial PZT materials [34|. A front matching layer is used to im- 
prove the sensitivity and bandwidth of the detector, and a back matching layer is used to 
170 improve the transmission of ultrasound into the backing layer and thus improve the band- 
width of the detector. The theoretical values for the suitable acoustic impedance of the 
matching layers (Zm) can be determined by 
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Zm — \/ ZiijZp (3) 

where Z^ and Zp are the acoustic impedance of water/tissue (or backing) and ferroelectric 
175 ceramic, respectively. The material properties of the piezoelectric material, the front match- 
ing layer, back matching layer and the backing are carefully chosen based on their acoustic 
impedances, and all important acoustic properties for these materials are listed in Table [H 

B. Functional and test models 

1. First functional model 

180 To have a rough assessment of the required characteristics of the detector such as center 
frequency (around 1 MHz) and bandwidth (>80%), the thicknesses of the active material. 



front and back matching layer and backing are determined using the KLM model |42|(See 
Sec. IV). The first functional model is then manufactured as shown in Figure [2](a), with 
dimensions of each layer listed in Table [Tll 



185 



2. Test models for minimizing radial resonances 



The large lateral dimensions of the detector result in radial resonances that interfere with 
the thickness resonance. To reduce the radial resonance of the element, the detector requires 
to be sub-diced into smaller units, which are acoustically isolated by air kerfs but electrically 
grouped by common electrodes to form a single composite element 43| (Figure [2](b)). To 

190 identify optimum lateral dimensions of the subdiced unit, five square- shaped, bare ceramics 
samples (CTS 3203HD) all having a thickness of 1.625 mm were manufactured as shown in 
Figure |2]^c). The different sizes represent the possible subdicing sizes applied to the 5 mm 
X 5 mm element to suppress the radial resonances. 3D finite-element method based models 
for these PZT samples are built using PZFlex (Weidlinger Associates Inc, Los Altos, CA) 

195 (See Sec. IV) to study the effect of sub-dicing, and to design the size of the sub-diced small 
unit. 



3. Second functional model 

This version results from the experiences with the first functional model and the test 
models above, and attempts to minimize interference from lateral resonances. For this, the 
200 first functional model is sub-diced into 0.9 mm x 0.9 mm units using a 100 /xm dicing saw. 
Only the front matching layer and active layer are sub-diced. The 25 units are acousti- 
cally isolated by air kerfs, but electrically grouped by two electrodes (Figure [2]^d)). The 
dimensions of each layer are listed in Table [ITl 



4- Final model 

205 With the second functional model above found suitable to minimize the presence of 
lateral resonance in the frequency range of interest, we studied the variation of passive layers 
dimensions on the frequency response using 3D FEM models. To increase the bandwidth 
of the detector, the thicknesses of the front and back matching layer need to be optimized. 
Three-dimensional FEM models are simulated to estimate the pulse-echo signals of the 

210 detectors, from which the frequency responses of the detectors are calculated to arrive at 
an optimized implementation. We manufactured a final model according to the geometrical 
parameters achieved as shown in Figure [2|^b,d). The dimensions of each layer are listed in 

8 



Table HI 



IV. NUMERICAL AND EXPERIMENTAL METHODS 



Simulation methods used 



1. ID KLM model 



An equivalent circuit based ID KLM model 



42| is used to obtain a rough assessment of 



the required detector performances such as center frequency and bandwidth. 



2. 3D FEM model 

220 For a more accurate estimation, coupled partial differential equations for piezoelectricity 
and acoustic wave propagation through passive material layers of the detector and coupling 
medium require to be solved 43|. This can be appropriately done by using 3D FEM models. 
In this study, 3D FEM models are built using PZFlex (version 3.0, Weidlinger Associates 
Inc, Los Altos, CA) as shown in Figure [2]^a) and (b). The finite element size is defined as 

225 1/30 of the ultrasound wavelength at the designed transducer center frequency (1 MHz). 
The roughly estimated results from KLM models are used as the starting point for the 
subsequent iterations between 3D simulations of detector performance and experimental 
characterization of functional models. The material properties used for FEM simulation are 
listed in Table [H 



230 B. Detector characterization methods 

1. Electrical impedance 

In order to ascertain the resonance characteristics, and to study the behaviour of the 
detector in an electrical measurement chain, the complex electrical impedance is measured 
using an impedance analyzer (4194A/B, Hewlett-Packard, Palo Alto, California). 

9 



235 2. Acoustic frequency response 



Two methods were used to measure the frequency response: (1) a transmit mode using 
a hydrophone, and (2) a pulse-echo mode. 

In method 1 (Figure [21^a)), the detector immersed in a deminerahzed water bath is driven 
by a broadband ultrasonic pulser /receiver (Panametrics 5077PR), and the generated pressure 
240 is probed in the far-field using a calibrated broadband needle hydrophone (0.2 mm diameter, 
bandwidth upto 10 MHz, Precision Acoustic Ltd., Dorchester, UK) with a known receiving 
transfer function Hq^/) . The frequency response of the detector (if(/)) is then calculated 
by: 



FFT{F(t)} 
^o(/) 



(4) 



Method 2 (Figure [3] (b)) uses the detector in pulse-echo mode, driven by the pulser- 
receiver. A stainless steel plate is placed in the far-field of the detector as an acoustic 
refiector, and the refiected signal (pulse-echo) is measured by the detector. According to 



44 



451], the frequency response of an 



the principle of reciprocity for a piezoelectric detector 
ultrasound transducer is the same when used as receiver or transmitter. Thus the frequency 
250 response of the detector is calculated by the square root of the FFT of the measured pulse- 
echo signal. 



3. Directivity 

The directivity of the detector is measured in transmit mode. The detector is mounted 
in a deminerahzed water bath and driven by the broadband pulser /receiver. The calibrated 
255 needle hydrophone is rotated in the far-field (60 mm) centered on the detector element to 
probe the emitted ultrasound field through 180°. For each scanning position, the peak- 
to-peak value is determined from the ultrasound pulse recorded by the hydrophone and 
plotted as a function of scanning angle, giving the directivity of the transducer at its center 
frequency. 

10 
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4- Sensitivity and minimum detectable pressure 



We measured these parameters using a modified substitution method [9] as outhned. We 



insonify the detector element using an ultrasound transmitter. The pressure incident on the 
detector element is progressively reduced, by reducing the voltage input to the transmitter. 
At each incident pressure, the voltage output of the detector is noted while averaging multiple 

265 times. This is continued till the lowest input possible. This detector output is defined 
as the signal. In a next step we note the noise voltage on the detector element without 
averaging and with no pressure incident on the element. The signal is plotted against 
pressure input and the trend is extended to intersect the voltage noise floor. This intersection 
point signifies SNR=1, and the pressure at which this signal is expected/obtained is the 

270 minimum detectable pressure. 

In the experiments we used a 1 MHz unfocused broadband transducer (V303, Panamet- 
rics) as a transmitter. Signals were processed by a prototype low noise pre-amplifier based 
on Analog Devices ADA4896-2 as the analog front-end and acquired using a high-speed digi- 
tizer NI-5752 (National Instruments). A calibrated needle hydrophone was used to ascertain 

275 the transfer function of the transmitter and convert the input voltage of the transmitter to 
known pressure that insonified the detector element. 

C. Imaging quality simulation 

To study the resolution and visibility of objects with a photoacoustic tomography (PAT) 
system employing the optimized detector, and comparing performances using detectors de- 
280 scribed in the literature, numerical simulations are performed using k-wave Matlab Tool- 



For the forward simulation, a 2D initial pressure distribution map (1024 x 1024 grid, 20 
cm X 20 cm size) is assigned in a 2D tomographic configuration. Three disc-shaped objects 
with diameters of 10 mm, 2 mm and 0.5 mm are located in the center region of the map. 
285 Homogenous initial pressure (value 1) is assigned to the objects, a pressure value of zero is 
given to the rest of the map. Homogenous acoustic properties are assigned to the medium 
(speed of sound: 1500 m/s, acoustic attenuation: 0, density: 1000 kg/m^.) Those pressures 
propagating outward are detected by a 5 mm detector rotating around the objects with 
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radius of 10 cm covering 360° with step size of 2". The time-domain photoacoustic signals 
290 are averaged over the surface of the detector at each detection position and saved for further 
processing. 

Three detectors are simulated, from: 

(1) The final model of this work. 

(2) The Kruger group (1 MHz center frequency and 100% bandwidth) 37l |. 

295 ( 3) The Oraevsky group (1.25 MHz center frequency and a bandwidth approaching 200%, 



derived from 



11 



19|). 



For (1), signals received by the detector are convolved with measured impulse response 
of the final model. For (2) and (3), the received signals are filtered using Gaussian bandpass 
filters with the corresponding center frequency and bandwidth. Correspondingly, three 
300 images are reconstructed by time-reversal of the processed signals. 



V. RESULTS 

A. First functional model performance 

Figure |4](a) shows the measured electrical impedance of the first functional model. A 
fundamental thickness resonance at 1.2 MHz together with a series of harmonic thickness 
305 resonances at higher frequencies, and a strong radial resonance at 330 kHz together with a 
second harmonic radial resonance at around 700 kHz, can be observed. The fundamental 
thickness resonance corresponds to the designed thickness of the active layer (1.625 mm) and 
the low frequency radial resonance arises due to the large lateral dimensions of the active 



layer (5 mm x 5 mm) 39 1. 



310 Figure IH^b) shows the measured frequency response of the first functional model using 
method 1 described above. Two peaks can be observed in the frequency domain response: a 
1.2 MHz peak caused by the thickness resonance, and a peak at 330 KHz caused by the radial 
resonance. Both peaks match with the measured electrical impedance peaks (Figure 111(a)). 
The first functional model has a center frequency of 1.2 MHz (maximum peak), with a -6 

315 dB bandwidth of 0.8 MHz. The fractional bandwidth is 67%. Due to the radial resonance, a 
secondary pulse (reverberation-like signal) after the main pulse is visible in the time domain 
signal, which is manifested as a passband ripple in the frequency domain. 

12 



The measured strong lateral resonance of the first functional model is not desirable. 



First the detector sensitivity is reduced by the lateral mode resonance 



43|. Second, strong 



320 lateral mode causes additional reverberation-like signals (Figure SJ^b)), which will adversely 



influence the image quality of the system 



46|. 



B. Test models performances 

Figure [5] shows the measured and simulated electrical impedance of the bare piezoce- 
ramic samples with thickness of 1.625 mm and lateral dimensions ranging from 5 mm x 5 

325 mm down to 0.5 mm x 0.5 mm. In general, FEM simulations exhibit good agreement with 
measurements (due to the manufacturing and measurement difficulties for the smaller bare 
ceramics, electrical impedance measurements for the 0.5 mm x 0.5 mm ceramic is not avail- 
able). The fundamental thickness resonance peak for all samples is always located around 
1 MHz, while the lateral resonance frequency peaks move towards higher frequencies as the 

330 lateral dimensions decrease. For lateral dimensions of 1 mm x 1 mm the radial resonance 
mode is suppressed. Further reduction in lateral dimensions to 0.5 mm x 0.5 mm does 
not significantly change the amplitude and location of the thickness resonance peak (Fig- 
ure [5](e) and (f)). Since more active material is lost with finer subdicing, leading to detection 
sensitivity loss, we prefer 1 mm x 1 mm as the final choice. 



335 C. Second functional model 

The measured electrical impedance of the second functional model in water is compared 
with 3D FEM simulation in Figure [6](a). Both simulation and measurement show that 
the lateral resonance is not visible in the electrical impedance curves as expected. The 
fundamental thickness resonance is located at around 1 MHz. 

340 The frequency response of the detector measured using method 2 compared with 3D FEM 
simulations is shown in Figure Et^c). Both measurement and simulation shows the detector 
has a center frequency of 0.9 MHz, however with a -6 dB fractional bandwidth of 48%, which 
is low for photoacoustic applications. 

13 



D. Final model 



The starting point for the thickness optimization is based on the 1-D KLM model j42|; this 
gives the optimized front matching layer thickness (^ml-f) of 0.58 mm and back matching 
layer thickness {Iml-b) of 0.54 mm as used for the first functional model (Figure 111(b)). 
In the first phase of the optimization, Iml-b is kept constant (0.54 mm), and Iml-f is 
tuned to obtain the largest bandwidth (Figure [71^a)-(c)). The optimum tML-F value is then 
350 determined (0.70 mm). The second phase is to keep the optimized tML-F constant, and to 
optimized tML-B- After the two phase optimization, the optimum Iml-f and tMi-B are 
determined (Figure [7](d)). 

The simulation results in Figure [7] show that the bandwidth of the detector is increased 
from around 50% (Figure [Tl^a)) to around 70% (Figure[7](c)) by optimizing the front matching 
355 layer. Further optimization of the back matching layer increases the bandwidth of the 
detector to more than 80% (Figure [7[d)), which is generally high for a PZT detector with 
single front matching layer [47 1. 



1. Frequency response 

The measured frequency response and time domain signal in the farfield (60 mm away) 
360 show good agreement with simulation results in Figures E^a) andlH^b). Compared to the 
ultrasound pulse for the first functional model (Figure Hl^b)), the optimized pulse is consid- 
erably shorter in length. The measured center frequency is 1 MHz, which is slightly higher 
than the model predicted 0.9 MHz (Figure |8]^b)). The measured -6 dB bandwidth ranges 
from 0.4 MHz to 1.25 MHz, which is also slightly larger than the simulation result (0.55 
365 MHz - 1.25 MHz), while the fractional bandwidths are similar for both (80%). 



2. Directivity 



Figure W,c) shows the measured directivity of the final model. The measurement results 
compare well with simulations, both showing a -6 dB acceptance angle around 20° at 1 MHz. 
When a 5 mm diameter hemispherical acoustic lens is placed on top of the detector, the 
370 directivity angle of the detector is expected to be enlarged to around 60° |31|. 
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3. Sensitivity 



The end-of-cable minimum detectable pressure of the detector-electronics was estimated 
to be 0.5 Pa. As seen the Figure [Hl^d)) the signal trend is extrapolated to intersect the noise 
floor to provide the point at which SNR=1. The corresponding pressure can be read off on 
375 the X-axis. 

E. Imaging quality 

The reconstructed images from the three cases are shown: for our detector in Figure [9]^b), 



Kruger system [37| in Figure IHI^c), and the Oraevsky system [ll|, ll9|] in Figure [Hl^d). It is 
clearly shown that Oraevsky detector, due to the superior bandwidth of the detector, pos- 
380 sesses excellent image quality and faithfully recovers the original initial pressure distribution. 
Ours and Kruger's system faithfully recover the 2 mm object. The 10 mm object is visu- 
alized with edge enhancement. The resolution of three systems is around 1-2 mm, which 
can be estimated from the reconstruction of the sub-resolution 0.5 mm object (see profiles 
in Figure El^e)). 

385 VI. DISCUSSION 

The final optimized single-element detector has a fractional bandwidth of 80%, which 
matches with the model prediction. The center frequency of the detector is around 1 MHz 
and -6 dB bandwidth is from 0.40 MHz to 1.25 MHz. From simulation of system resolution 
and image quality (Figure IHl^b)), it is found that the object with diameter of 2 mm can 

390 be faithfully recovered, and the 10 mm object can be visualized with acceptable distortions 
(edge enhancement and ring-shaped artifacts). The resolution of our system reaches the 
designed goal of 1-2 mm (Figure |9]^e)). The simulations are performed with a high optical 
contrast of the object, which gives large signal-to-noise ratios. This condition is applicable 
for a sensitive detector like ours: frequencies beyond the -6 dB bandwidth of the detector 

395 will still contribute to signal detection, which makes the resolution of our system better 



than the value calculated from Eq.(2) 19|] taking only -6 dB bandwidth of the detector 
into account. The resolution metrics used via k-wave simulations do not necessarily line 
up with those from Eq.(2) because Eq.(2) is derived only from the frequency contents of 

15 



the photoacoustic signals generated from objects with varying sizes. This approach gives 
an indication for designing the ultrasound detectors, while imaging approach via k-wave 
simulation is suitable to characterize the overall system resolution. This indicates that the 
present detector is suitable to detect tumors during the beginning of the vascular phase. 
The image quality can be further improved by enlarging the bandwidth of the detector 



using two or more front matching layers 



However, a larger bandwidth reduces the 



405 sensitivity of the detector, and increases complexity and expense during transducer design 
and development and in final production. 

The directivity angle of the transducer is around 20°, which can be increased to around 
60° using an acoustic lens. For simplicity, considering the breast as a hemisphere with 
diameter of 10 cm, the transducer with directivity angle of 60° is required to be placed only 

410 5 cm away from the breast. This distance is suitable for the application. 

We have consolidated the most important output characteristics of the optimized ultra- 
sound detector in Table llllt in row 7 described as the PAM-II system. We have also provided 
specifications of various ultrasound detectors described in literature for comparison. The 
MDP value of our detector is 0.5 Pa, which is 160 times lower than for the detector used 

415 in the first generation of the Twente Photoacoustic Mammoscope previously developed in 



our group [20| . To the best of our knowledge, this is the lowest measured MDP reported for 
a detector in photoacoustic breast imaging (Table [In|) . The total noise performance of the 
system could possibly further be improved with better electrical shielding and grounding. 

The bandwidth of our detector is 80%, which is slightly broader than the detectors used 
by Kruger group and Kyoto University in their breast imagers (See Table UTTj) . The detector 
from the Oraevsky group has an ultrabroad bandwidth approaching 200% (estimated from 



Refs. 
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19|). which provides impressive image quality as shown in Figure [9]^d). The image 



quality of our detector can be improved with the use of a deconvolution operation as is 



performed by several groups |49|-|52| to compensate for finite bandwidths effects. Further 
425 breast tumors are known to be heterogeneous with a scattered distribution of absorbing 
regions. In such a case our detector is eminently suited to faithfully image a collection 
of small absorbing structures, which makes up the tumor mass. We are aware that there 
is room for improvement, nevertheless our design strategy has resulted in a detector with 
acceptable bandwidth which is well suited for sensitive clinical breast imaging due to its 
430 unprecedented high sensitivity. 
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Detector arrays will be manufactured based on the single-element detector described in 
this work for use in the second version of the Twente Photoacoustic Mammoscope (PAM- 

11) bi- 



VII. CONCLUSION 

435 A single-element PZT, large-aperture, sensitive and broadband detector is designed and 
developed for photoacoustic tomography of the breast. Finite-element based models are used 
to optimize an initial detector to reduce the radial resonance and optimize the bandwidth of 
the detector. The center frequency and -6 dB fractional bandwidth of the optimized detector 
are 1 MHz and around 80%, respectively. The minimum detectable pressure is 0.5 Pa, which 

440 is more than two orders of magnitude lower than in our first generation photoacoustic breast 
imaging system and among the lowest reported in the literature. Detector arrays will be 
manufactured based on the design of this single-element detector for use in the second version 
of the Twente Photoacoustic Mammoscope. 
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LIST OF FIGURE CAPTIONS 



Figure 1: Schematics of a 2D photoacoustic tomography system showing the necessity for 
the acceptance angle to encompass the object for coherent signals detection from all angular 
position in performing reconstruction. 

Figure 2: (a) Schematic of first functional model, (b) Schematic of subdiced second func- 
595 tional model and final model. (c) Photograph of the bare PZT samples with different lateral 
dimensions (label below). Two triangular samples are shown in the photograph, however, no 
related result is reported in this work, (d) Photograph of the final single-element model. 

Figure 3: Schematics of the setups for detector frequency response measurements, (a) 
Transmit mode, (b) Pulse-echo mode. 
600 Figure 4: First functional model performance, (a) Electrical impedance of the first functional 
model measured with water load, (b) Measured transmission impulse response and frequency 
transfer function of the first functional model using a hydrophone, (time domain left axis, 
frequency domain right axis) 

Figure 5: Test samples of PZT: measured and simulated electrical impedance in air with 
605 lateral dimensions (a) 5 mm x 5 mm; (b) 4 mm x 4 mm; (c) 3 mm x 3 mm; (d) 2 mm x 2 mm; 
(e) 1 mm x 1 mm and (f) 0.5 mm x 0.5 mm. No measured impedance available for 0.5 mm x 
0.5 mm PZT due to the practical limitations in manufacturing and measuring. 

Figure 6: Second functional model performance: (a) Measured and simulated electrical 
impedance in water, (b) Measured and simulated pulse-echo ultrasound signal. The reflector is 
610 placed in the far-field of the detector. The time delay is removed, (c) Measured and simulated 
frequency response. 

Figure 7: Towards optimized final model: simulated pulse-echo signal and the frequency 
response of the sub-diced detector with different front- and back- matching layer thicknesses, 
(a) tML-F- 0.58 mm, Iml-b'- 0.54 mm. (b) tuL-F'- 0.55 mm, Iml-b'- 0.54 mm. (c) Iml-f'- 
615 0.70 mm, Iml-b'- 0.54 mm. (d) tML-F'- 0.70 mm, Iml-b'- 0.48 mm. 

Figure 8: Final model performance: (a) Measured and simulated far-field pulse of the final 
model time-shifted to origin. For the measurement, the pulse is probed using a calibrated broad- 
band needle hydrophone in the far-field at distance 60 mm, on center axis, (b) Measured and 
simulated frequency response of the detector, (c) Measured and simulated directional sensitivity. 
620 (d) Measured sensitivity and minimum detectable pressure (MDP). 
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Figure 9: (a) Initial pressure distribution used in the forward simulation. The gray dashed 
circle indicates the detector scanning positions, (b) Reconstructed image using signals detected 
by the final model, (c) by the transducer with 1 MHz center frequency and 100% fractional 
bandwidth (Kruger et al 1999 37j) and (d) by the transducer with 1.25 MHz center frequency 
625 and 200% fractional bandwidth (Andreev et al 2003 [19]and Ermilov et al 2009 (nl). (e)Profiles 
at position X=0 mm from the initial pressure distribution in (a) and reconstructed images from 
(b),(c) and (d). 
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LIST OF TABLE CAPTIONS 



Table 1: Properties of the materials used for the detector in 3D FEM simulations. Properties 
630 of the PZT material are from reference 0, 41|, and the properties of the matching and backing 
layers are from reference j4l| . 

Table 2: Layer thicknesses of functional models (Figure [21(a) and (b)), each layer has a 5 
mm X 5 mm square-shape surface. 

Table 3: List of US detectors used by different groups in the photoacoustic (thermoacoustic) 
635 systems for breast imaging. 
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TABLES 



TABLE I. 



Properties Matching layers Active layer Backing 

Material Electrical conductive epoxy PZT(CTS 3203HD) Elastosil 

Impedance (MRayl) 6.5 38 3.4 

Density (kg m-^) 3140 7800 1870 

Longitudinal velocity (m s^-*^) 2068 - 1818 

Shear velocity (m s~^) 994 - 873 

Elastic compliance - sf^=1.56, s^=-0.420, - 

(m^/NxlO-^i) sf3=-0.823, sf3=1.89, 

Piezoelectric strain coefficient - di3=2.95, ^33=6. 64, 

(mVNxlO-i°) di5=5.60 

Relative permitivity - Kfi=2417, iC|^=3331 - 

Dielectric loss - 0.028 

Mechanical quality factor - 66 - 
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TABLE 11. 

Layer description Material 1®* and 2"^^^ functional models Final model 



thickness (mm) thickness (mm) 

Front matching layer Electrical conductive epoxy 0.590 0.700 

Active layer CTS 3203HD 1.625 1.625 

Back matching layer Electrical conductive epoxy 0.788 0.480 

Backing Elastosil 10 10 
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TABLE III. 



VJi KJ U.]JO 




Detector Element 


Active 


Center 


U V V 


MDP 






1 I 
elements 


geometry material 


frequency 




[fa) 








t lllill 1 




(MHz) 






(1) Kruger group 


TA [17] 




13 


1-3 composite PZT* 1 


/ U/O 




(2) Kruger group 


PA [lOj 


128 


3 


1-3 composite PZT* 5 


70%* 




(3) Wang group 


TA&PA [12] 


1 


13/6 


Piezocomposite^^ 


2.25 


60-120%t 




(4) Oraevsky group 


LOIS-64 [IJJ 


64 


20 X 3 


PVDF 


1.25 


170% § 


** 


(5) Kyoto University 


PAM [48] 


345 


2x2 


Piezocomposite 


1 


>70% 




(6) University of Twente PAM [24] 


590 


2x2 


PVDF 


1 


130% 


80 


(7) University of Twente PAM-II^ 


1 


5x5 


CTS 3203HD 


1 


80% 


0.5 



* R. Kruger (2012), private communication. 

^ L. V. Wang and M. Pramanik (2012), private communication. 

^ A. A. Oraevsky (2012), private communication. 

** Various values are reported in literature including measured [2]|, estimated |l9l . 
measured MDP of 1.8 Pa from A. A. Oraevsky (2012), private communication. 
Planned. 
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